Interior Gap Superfluidity 
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We propose a new state of matter in which the pairing interactions carve out a gap within the 
interior of a large Fermi ball, while the exterior surface remains gapless. This defines a system 
which contains both a superfluid and a normal Fermi liquid simultaneously, with both gapped and 
gapless quasiparticle excitations. This state can be realized at weak coupling. We predict that a cold 
mixture of two species of fermionic atoms with different mass will exhibit this state. For electrons in 
appropriate solids, it would define a material that is simultaneously superconducting and metallic. 



Recent developments in ultracold alkali atomic gases 
[jlj have revitalized interest in some basic qualitative 
questions of quantum many-body theory, because they 
promise to make a wide variety of conceptually inter- 
esting parameter regimes, which might previously have 
seemed academic or excessively special, experimentally 
accessible. With this motivation, and stimulated by ques- 
tions in quantum chromodynamics (QCD) at high den- 
sity ||, ||, U , we here revisit the question of fermion pair- 
ing between species whose Fermi surfaces do not precisely 
match. We have found a possibility that seems to be new 
and certainly is interesting, and which could turn out to 
be relevant even for conventional solids. 

The standard Bardeen-Cooper-Schrieffer (BCS) || 
theory of superconductivity describes pairing between 
particles of equal and opposite momentum near a com- 
mon Fermi surface. For classic s-wave superconductors 
the pairing occurs between electrons of opposite spin. In 
the presence of a weak magnetic field, and in particular 
in the case of ferromagnetic order, the Fermi surfaces of 
the opposite spins will not match, and the Cooper pairing 
instability, which was enhanced by vanishing energy de- 
nominators, will no longer occur at arbitrarily weak cou- 
pling. Larkin and Ovchinnikov and independently Fulde 
and Ferrell Q showed that in this circumstance it might 
be favorable to effectively relatively translate the Fermi 
surfaces, pairing at a non-zero total momentum (LOFF 
phase) . 

A simpler situation, conceptually, is that pairing oc- 
curs between two species whose Fermi surfaces do not 
match simply because their densities or effective masses 
differ. This possibility arises in several contexts, (i) In ul- 
tracold atom systems, it could occur simply because there 
are atoms of different elements, (ii) In solids it could oc- 
cur for electron populations in two different bands, (iii) 
In QCD it occurs for different species of quarks (up, 
down, strange). If the mismatch is small and the two 
species are alternative states of the same particle (such as 
the spin up and down states of electrons or two hyperfme- 
spin states of cold 40 K or 6 Li atoms as prepared in ex- 
periments [0, ||, ^|, |l(|), it can be favorable to equalize 
the Fermi surfaces, absorbing a cost in kinetic energy, 
and then to pair at zero momentum following BCS. For 



larger mismatches, LOFF-type ordering can occur. More 
elaborate forms of position-dependent ordering, with the 
superfluid gap having standing wave or even crystalline 
structure, have been found to be favorable in models of 
QCD at high density Q. 

None of these possibilities, however, extrapolates to 
what we might expect at strong coupling. Given strong 
attraction between the species, we would just expect to 
bind as many paired quasi-molecules as possible. At 
low temperature, these quasi-molecules will form a Bose- 
Einstein condensate. The residual unpaired particles will 
constitute a separate normal fluid. It is natural to inquire 
whether there is a weak-coupling phase that matches this 
behavior qualitatively. We now identify such a phase. 

Consider a homogeneous fermion gas in three dimen- 
sions containing two species (a = A, B) obeying simple 
parabolic dispersion relations, described by the Hamilto- 



pa 



pp'q 



-0t ^ 



A, 



q+p' 



(1) 



where e a (p) = p 2 /2m Q — fi a . See Fig. |]. (More precisely, 
we assume that this interaction exists so long as the mo- 
menta of the particles are within a strip of size A around 
the smaller Fermi surface; A will later serve as an ultra- 
violet cutoff.) Our heuristic analysis will not distinguish 
whether or not the species are strictly conserved sepa- 
rately. We define chemical potentials so that the Fermi 
surfaces for both species are at tp = 0. We shall be 
interested in cases where wia < tub and > Mb, in 
such a way that the Fermi momentum for the species B 
is greater than that for the species A, pp > Pp. 

We suppose that there is an attractive effective inter- 
action in s-wave between particles of different species 
(g < 0), and that the coupling is weak, so that we 
can construct our ground-state by modifying the ground 
state of the non-interacting system. If the Fermi surfaces 
matched, the attractive interaction would trigger stan- 
dard BCS superfluidity, with Cooper pairing of equal 
and opposite momenta. The BCS wave function, how- 
ever, postulates either zero or double occupancy of the 
paired modes, and it is incompatible with keeping the 
modes of species B between pp and pp completely filled. 
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FIG. 1: The prototype situation where we anticipate forma- 
tion of an interior gap superfluid at weak coupling. There are 
two species of fermions with different band structures, here 
both taken as isotropic and parabolic, but with different ef- 
fective masses and different sizes in momentum space. When 
the larger Fermi ball derives from relatively flat band (large 
effective mass) and the interaction near the momentum sur- 
face defined by the smaller Fermi sphere is attractive, it can 
be favorable to form correlated pairs near this smaller sphere, 
even at the cost of promoting some particles of the heavier 
species to the exterior Fermi sphere. One will then have both 
superfluidity, with a momentum gap at the smaller sphere, 
and normal Fermi liquid excitations at the larger sphere. 



If we are to support pairing of total momentum zero we 
must promote some particles of species B up to momenta 
near pf , thus carving an interior "trench" of the species 
B Fermi sea near momentum p = pp . 

There is competition between the energetic cost of such 
promotions and the gain from pair-formation, and it may 
not be obvious whether there can be a net profit in any 
non-trivial case. To assess this, let us suppose that the 
pairing introduces a momentum gap of order k. By this 
we mean that in an interval of order n around pp we will 
take superpositions of unoccupied and doubly occupied 
states, as in ordinary BCS theory. In particular, we do 
not automatically fill the single-particle states for species 
B, even though they are below the free-particle Fermi 
momentum pf. One could also, more awkwardly but 
perhaps more properly, speak of normalized energy gaps 
of order npf /m a for the two species. However phrased, 
the point is that it is important to prepare an equal num- 
ber of modes to pair, and state-counting takes place in 
momentum space. The condensation energy must be of 
the same energy as the spectral displacement of the par- 
ticles, so we have for the energy gain e pa i r per pair e pa i r ~ 
Ppn/fh, with ?7i = mArriB / (m,A + mB) the reduced mass. 
On the other hand the density of pairs n pa i r is of order 
"■pair ~ 'iirpp 2 n. In order to accommodate the depletion 
of species B, which is of the same order as the number of 
pairs, we must promote a corresponding number of parti- 
cles from pp to pp, which costs 2 ^ (pf —pp ) per pair. 



Putting it all together, to make a net profit we require 
Epair ~ PpK/fh > (pf - pj, )/2m B . 

The region of pairing in momentum space will be 
strictly in the interior of the larger Fermi surface - 
clearly distinct from its boundary — and we will real- 
ize interior gap superfluidity, if pf—Pp > k. This is 
compatible with our earlier condition for 

1> P f+P$ m A 

2pp mA + rriB 

This consistency condition can be satisfied, specifically 
for ttib ^> mA- Note particularly that it is independent 
of the gap k. Looking back to our net profit condition, 
we see that k can be taken arbitrarily small, for suf- 
ficiently large ttib/to^. Thus interior gap superfluidity 
can take place at weak coupling, where the mean-field as- 
sumptions implicit in this heuristic analysis are valid. In 
the class of models under discussion, therefore, we find 
a robust weak-coupling phase characterized by a (mo- 
mentum) gapped Fermi surface interior to a surface with 
unpaired excitations. For charged fermions, it is a phase 
that is simultaneously superconducting and metallic at 
zero temperature. 

To construct the interior gap ground state explicitly 
we generalize the standard BCS wavefunction as follows 

|* IG )= 1] (sin0 p +co8Ml p 4,-p) n ^ P i°) 

P|<PA |p|>PA 

where the p 's and p/\ are variational parameters. As 
usual, there is a manifold of degenerate states fea- 
turing an overall relative phase between the sin# p 
and cos6* p factors. The order parameter is of usual 
form: (ipAp^B - p )ig = sin# p cos0 p . Upon varia- 
tion with respect to 9 P and pa , we find cos 2 9 p = 

K^vtfM and A = ^ f + pf) - 

ihl6A 2 m A m s + (pf 2 - pf f} 1 / 2 , with e+ = §( e £ + 
£p). The gap parameter, defined here as A = 
~9 E| P |< PA (V'Ap' ! /'L - p )ig, satisfies the integral equation 
1 = -s£|pI<pa^T^- 

An important departure from standard BCS theory oc- 
curs because the energy difference between paired and 
unpaired modes no longer becomes arbitrarily small. For 
this reason one is not doing degenerate perturbation the- 
ory, and does not encounter a true infrared divergence 
(vanishing energy denominators) in the Cooper pairing 
channel. As a consequence, the gap equation supports a 
non-zero solution only for \g\ > g c , with 

9c t \ — -y , (3) 

where we have introduced the generalized density of 
states N+(Q) = J2 P ^( € p) ancl Po is the point where 
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e+ (Po) = 0. Note however that g c — > when m^/m^ — > 
oo for fixed p F ' , so that interior gap superfluidity can 
be favorable for arbitrarily weak attractive interactions, 
as we anticipated. Numerically, we typically find that g c 
is rather smaller. 

It is straightforward to calculate the condensation en- 
ergy. Up to terms of higher order in A, we find 



Ei G -E N = -i7V+(0)A 



XAe 



arcsinh x / 



(4) 



where x\ = 6 + (pa)/2A and is the normal state en- 
ergy at A = 0. For weak pairing pa is close to pp, in 
which case £ + (pa) is negative and finite. In the limit 
A — ► 0, ia approaches minus infinity, and the two terms 
in the brackets cancel. One sees that, for a given gap pa- 
rameter A, the interior gap state gains less condensation 
energy than a conventional BCS state. 

Thus we have demonstrated that the normal state 
is unstable against formation of an interior gap super- 
fluid when the attractive coupling is strong enough, i.e., 
\g\ > g c . The LOFF state is another candidate for pairing 
of mismatched Fermi surfaces. We have calculated the 
ground-state energy difference between these two candi- 
date states numerically. Fig. || shows a typical phase dia- 
gram. We have considered here only the simplest LOFF 
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FIG. 2: A typical phase diagram, with the mass ratio 
tub/hia = 7. Our coupling constant is g = —iivh 2 a s /m, 
where a s is the scattering length for s-wave scattering be- 
tween species A and B, and n t = n A +n B is the total density 
of particles, calculated using the formula = (67r 2 n a ) _1//3 as 



kiPF+PF) 



for a free Fermi gas. ApF = Pf — Pf an d p^' 
The dash-dotted line is plotted using the approximate ana- 
lytical formula ([]), while the solid line is plotted by adapting 
Eq. (3.12) of Takada and Izuyama |L2j to the case at hand. 
Were it plotted directly from our numerical simulations, the 
phase transition line between interior gap and LOFF would 
be slightly lowered. 

state, with pairing order parameter oc Aloffc 1< ^ x with 
|Q| of order pf-p F § (T|. 

A few remarks concerning this phase diagram are in or- 
der. The critical coupling constant g c in Eq. (^|) in princi- 
ple governs the (second-order) phase transition from the 
normal state to the interior gap state, rather than that 



from the LOFF state. However, we find in our numer- 
ical calculation that the phase transition line between 
the interior gap and LOFF states roughly coincides with 
the onset of the interior gap ordering. Indeed, since an 
LOFF particle pair carries a total momentum Q ^ 0, the 
pairing process in the LOFF state occurs mostly within 
the intersection of two closed shells of thickness k, one 
centered on Q and the other on zero. The intersection 
of these shells is a closed ring of thickness k. By con- 
trast, for the interior gap state pairing occurs within a 
full two-dimensional shell. Since the density of states in- 
volved in pairing for the interior gap state is larger than 
that for the LOFF state, the former develops an order 
parameter which increases (as a function of intrinsic cou- 
pling strength) exponentially faster than the latter, and 
rapidly dominates once it sets in. On the other hand, the 
LOFF phase can be realized also for tua > tub, that is 
when the heavier fermion has the smaller Fermi surface, 
when the interior gap phase is not available. This is the 
case of primary interest for possible phases of quark mat- 
ter in neutron star interiors, where the B species is the 
strange quarks. BCS states correspond to the region in 
our phase diagram where the two species have approxi- 
mately equal Fermi momentum, ApF ~ 0. 

The average particle occupation number n^' B has an 
unconventional form, rip = rip = cos 2 9 p for |p| < p<\, 
and rip = and rap* — 1 for pa < |p| < Pf, where pf is 
the shifted Fermi momentum of species B particle due to 
pairing interaction: pp ~ pf (1 + pofhAe arcsinhx& /pp 3 ). 
Fig. ^| shows the average particle occupation number as a 
function of momentum at zero temperature. In the inte- 
rior gap state, the pairing correlation smears the species 
A Fermi surface slightly — a fraction of single-particle 
states are depleted below pp and inserted back between 
the normal state Fermi surface at pf, and the maximum 
pairing surface at pa- This distribution does not differ 
qualitatively from what one finds in a conventional BCS 
superconducting state. Species B displays a more dra- 
matic contrast. Some modes that are occupied in the 
normal state, in the Fermi ball interior around p = pp, 
are now depleted, and the deficit is made up at the top of 
the Fermi surface. This enlarges the Fermi surface from 
Pp to pp. Discontinuities in the distribution rip occurs 
at both pa and pp . 

Approximate quasi-particle excitations of the interior 
gap state are obtained by diagonalizing the Hamilto- 
nian (Q) at mean-field level using the Bogoliubov trans- 
formation 7^ p = sinflp^ — cos (JpV's.-p and 72 _ p = 
sin 6 p tp\ 



cosOptpAp- These operators create quasi- 
particle excitations 7J 2-pI^ig) from the interior gap su- 
perfluid state. Their spectra are given by 



£i, 2 (p) = 50£- e £)±V e ? +A2 > ( 5 ) 

with all energies measured from the shifted Fermi sur- 
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FIG. 3: a and b: A plot of average particle occupation num- 
ber rip and rip, respectively, in the interior gap superfluid 
state in comparison with those in a conventional BCS state. 
For the comparison, we assume the same gap parameter and 
a matched Fermi surface at p = pp. c. Quasiparticle en- 
ergy spectra to add a species A or B particle to the system. 
Corresponding quasi-particle states are V'apI^'ig) (' — ' line) 

and V'bpI^ig) (' ' line), d. Quasiparticle energy spectra 

to remove a species A or B particle from the system. Cor- 
responding quasi-particle states are V'ApI^ig) (' — ' line) and 
^SpI^ig) (' ' line). 



faces defined by -pp ~ pp and pp . Unlike in conventional 
BCS theory, there are two branches of excitations: E\ (p) 
is gapped while E 2 (p) is gapless at both |p| = pa and 
|p| = pg. Notably, '— i?2 P ' becomes negative for mo- 
menta |p| between pa and pp. This is interpreted con- 
sistently as indicating that all such states are filled by 
species-B particles (note sin# p = 0), and that the abso- 
lute value |i?2p| is the energy to create a hole excitation 
there (Fig. |). 

Given the explicit form of the quasiparticles and their 
spectrum, phenomenological consequences can be derived 
along standard lines. The novelty of the interior gap state 
is that a large manifold of low-energy "normal state" 
excitations coexists with superfluidity. This spectrum 
could be probed directly in tunneling experiments. At 
finite temperature the normal state excitations will be 
excited, and the appropriate description will involve a 
two-fluid model incorporating dissipation. In these re- 
gards there is some resemblance to conventional superflu- 
ids whose order parameter has nodes, such as the p-wave 
superfluid in liquid 3 He, or the d-wave superconducting 
cuprates. But these states differ from interior gap super- 
fluids both quantitatively, in that the density of gapless 
modes is much smaller, and qualitatively, in that they in- 
volve breaking of rotational symmetry Another partial 
analogue is the Abrikosov-Gor'kov gapless superconduc- 
tivity with magnetic impurities Jl3| ; but of course here 
we do have a gap, and impurities are not a central issue. 



Interior gap superfluidity will be realized in a two- 
species mixture of fermionic cold atoms with different 
mass. Recent theoretical Q [Te| [l?], |l|, |l9) and 
experimental ^ ||, [| [To| efforts point to the possibil- 
ity of superfluidity in two-state mixtures of 6 Li or 40 K 
atoms. A stable mixture with different mass could be 
realized, for instance, in the 6 Li and 40 K atomic gas 
p0| . Despite its qualitative difference from BCS-type 
states, the interior gap state is estimated to have a su- 
perfluid transition temperature of the same order as that 
of Refs. @ £l|, [16], [17], [l8[ [l9| , for not too weak coupling. 

Also, we perceive no problem of principle forbid- 
ding the realization of interior gap superfluidity in elec- 
tron gases, where the species are electrons from differ- 
ent bands, which can have markedly different effective 
masses. The case of electrons coming from two bands 
differs however from the atomic case in that one should 
specify a single density and an energy offset between 
band minima, instead of two independent densities (or 
chemical potentials). There can be phase transitions as 
a function of the single density, for example between in- 
terior gap and conventional BCS-type order. Also, of 
course, the dispersion relations can be different from sim- 
ple parabolas, which has interesting consequences. We 
shall return to these questions in a future publication. 

We have benefited from discussions with J. Bowers, C. 
Honerkamp, W. Ketterle, K. Rajagopal, and X.-G. Wen. 
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U.S. Department of Energy (D.O.E.) under cooperative 
research agreement #DF-FC02-94ER40818. 
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